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Abstract- The paper concentrates on the study of propagation 
of thermoelastic waves in a homogeneous, transversely 
isotropic, thermally conducting elastic plate bordered with 
layers (or half-spaces) of viscous liquid on both sides in the 
context of non classical theories of thermoelasticity. Complex 
secular equations for symmetric and antisymmetric wave 
motion of the circular plate, in completely separate terms, 
are derived. Finally, in order to illustrate the analytical results, 
the numerical solution is carried out for transversely isotropic 
plate of cobalt material bordered with water by using the 
functional iteration method. 

Index Terms- Crested waves, viscous, thermal relaxation, Biot's 
constant, transversely isotropic. 

I. Introduction 

The use of sandwich structures in various applications, 
especially in many engineering phenomena, including the 
response of soils, geological materials and composites, the 
assumption of isotropic behaviour may not capture some 
significant features of the continuum response. The 
formulation and solution of anisotropic problems is far more 
difficult and cumbersome than its isotropic counterpart. There 
are reasonable grounds for assuming anisotropy in the 
continents. Anisotropy in the earth crust's and upper mantle 
affects the wave characteristics considerably. Also, 
Ultrasonic Lamb waves have been investigated extensively 
for damage detection in advanced composite materials. 
Banerjee and Pao [3] investigated the propagation of plane 
harmonic waves in infinitely extended anisotropic solids 
taking into account the thermal relaxation time. 

Dhaliwal and Sherief [4] extended the generalized 
thermoelasticity theories to anisotropic elastic bodies. Zhu 
and Wu [5] brought out detailed analysis of the characteristics 
of Lamb waves in elastokinetics. Sharma and Singh [6], 
Sharma and Pathania [7] and Sharma and Walia [8] 
investigated the propagation of Lamb waves in thermoelastic 
solid plates bordered by liquid media maintained at uniform 
temperature. Recently, Sharma and Sharma [9] investigated 
the propagation of Lamb waves in thermoviscoelastic plate 
loaded with viscous fluid layer of varying temperature. 
Hussain, Ahmad and Asawa [10] investigated analytically 
and through experimentation the discharge capacity of a 
sharp-crested circular orifice in an open channel under free 
flow conditions. The present investigation is concerned with 
the study of circular crested wave propagation in an infinite 
homogeneous, transversely isotropic, thermoelastic plate 

©2012 AMAE 
DOI:01.DPIE.02.01.45 



bordered with an viscous liquid layers or half spaces on both 
sides.lt is noticed that the motion for circular crested waves 
is also governed by Rayleigh-Lamb type secular equations 
as in the case of a rectangular plate. More general dispersion 
equations of Lamb type waves are derived and discussed in 
coupled and uncoupled theory of thermoelasticity. Numerical 
solution of the dispersion equations for cobalt material has 
been carried out and presented graphically. 

H. FORMULATION OF THE PROBLEM 

We consider an infinite homogeneous, transversely 
isotropic, thermally conducting elastic plate of 

thickness 2d initially at temperature T Q , having a viscous 

liquid layer of thickness ]% on both sides. We take origin of 

the co-ordinate system (r, 0, z) on the middle surface of the 

plate and % -axis normal to it along the thickness. The plate 
is axi-symmetric with respect to £ -axis as axis of symmetry. 
We take r — z P lane as the plane of incidence. From Sharma 
and Pathania [7] and Chandrasekhar [1], we have taken the 
basic governing equations in non-dimensional form for 
generalized thermoelasticity in the absence of body forces 
and heat sources and for viscous liquid (homogeneous liquid) 
medium respectively. 
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Here t Q and t, are the thermal relaxation times. 8 ik is the 

Kronecker's with fc — \ for Lord-Shulman (LS) theory and 

k — 2 for Green-Lindsay(GL) theory of generalized 
thermoelasticity. s is the thermomechanical coupling 

constant and v, is the longitudinal wave velocity in the 

thermoelastic half-space. -/ t \ / n \ is the 

r u\r,z,t)=\u,\j,w) 

displacement vector, T (r , z , t ) is the temperature change, 

Tq is the reference temperature of solid, c, is the velocity 

of sound in the liquid, a * is the coefficient of volume thermal 

expansion, s, is the thermomechanical coupling in the liquid. 
Along the z direction for the top liquid layer ( j = l ) and 

for the bottom liquid layer () = 2 ) . C, , p , are velocity of 
sound, the density of the liquid respectively, X. is the bulk 
modulus, C v is the specific heat at constant volume of fluid. 
Tj is the temperature deviation of liquid medium from 
ambient temperature T 

u L (r, z,t)= [u L ,0, w L j ( j = 1,2) is the velocity 

vector. U, and W L are respectively. Along the direction for 

the top liquid layer and for the bottom liquid layer .Moreover 
liquid is assumed to be thermally non-conducting 
hypothetical solid. Here the dot notation is used for time 
differentiation and comma denotes spatial derivatives. 
In the liquid layers, in addition to introducing non-dimensionl 
quantities, we take 
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where <P L and \f/ L ,{j =1,2) are respectively, the scalar 
and vector velocity potential function, we obtain, 
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III. Boundary Conditions 

Boundary conditions at the solid-liquid interfaces i — +d 
are: 

r = + H (T-T L ) = (10) 
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where // is Biot's heat transfer coefficient 

IV. SOLUTION OF THE PROBLEM 
We assume solutions of the form 
{ u : m% T}={J, l> I VJ, fc J WJ, iir ) JSe* *--» <11 1) 

K JTi, j= {* tj '« (H^ '. fe*0fc* f —' ° (1 1.2) 

_ a 
where c - ~y is the non-dimensional phase velocity, CO is 

the frequency and E, is the wave number. 7 , /( are 

respectively, the Bessel functions of order zero and one. e 
is the angle of inclination of wave normal with axis of 
symmetry, m is still unknown parameter. V , W respectively 
are, the amplitude ratios of displacement w , temperature j 

and <J> L , *F L are the amplitude ratios of velocity potentials 

L ,w l respectively. 

Upon using solution (11.1) in (1) - (3). We obtain 

m 6 + Am 4 +Bm 2 +C = (12) 

Here 
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Equation (12) being cubic in m admits six solutions for 

m, with property m 2n = —fn 2n _ 1 , n = 1,2,3. On seeking 

the solution of (l)-(3) in the form (11.1), displacements, 
temperature change can be expressed as 

w " T.l-, v <(- A sin <f M ,< + S, ok fB.jtK^rV-*- 

where m are the roots of (12). A , fi are amplitudes. 
V ? , W are amplitude ratios which can be expressed as 
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where 

On using stress strain relation 

\ r J 

axial stress and shear stress in T Z plane is obtained as 

where 



'4^y*±. *=\ 3 fe+^K.ff-i.3 



•; — — ,*=5 



«!-*. 



(17) 



"i. 



,f=5 



On seeking solution of form (14.2) for (4) in liquid layers 
d <z<d + h and ~(d + h)< z<-d , such that the 

acoustical pressure vanishes at z — i\d + Jl), we obtain 
(*<, Tx, H 1 ^^^** !?-(<*+*)]«" (IS) 

T t , =A,smim J [z-(d + A)]e"" (20) 

^ L _ = - A n sin d= js 9 [i + (<i + A)]e "'" (21) 

where A t ,q = 7 , 8 , 9 , 10 are amplitudes. 



S. 



ftrf^d+f[)--to»l] 



"~^ 1 



(22) 



V. DERIVATION OF THE SECULAR EQUATIONS 

By invoking the interface conditions (13), we obtain a 
system of ten simultaneous linear equations. System of 
equations have a non-trivial solution if the determinant of 

the coefficients of A vanishes, which leads to a 

characteristic equation for the propagation of modified 
guided thermoelastic waves in the plate. The characteristic 
equation for the thermoelastic plate waves in this case, after 
applying lengthy algebraic reductions of the determinant 
leads to the following secular equations 
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Here 

T. = tan wi. 
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G 3 ,G 5 , G 7 and G 3 ,G 5 , G 7 can be written from 
G, andG,' by replacing subscript (3,5,7) of elements in 

G x and Gj' with (l , 5 , 7 ), (l , 3 , 7 ), (l , 3 , 5 ) respectively with 

the exception that (l,3), (2,3) element in both 

G 7 and G 7 are with positive sign and (3,3) element is given 

by W s [m s - HT 5 ) .The superscript +1 corresponds to 
asymmetric and -1 refers to symmetric modes. 

The Rayleigh-Lamb type equation also governs circular 
crested thermoelastic waves in a plate bordered with layers 
of viscous liquid or half space viscous liquid on both sides. 
Although the frequency wave number relationship holds 
whether the waves are straight or circularly crested, the 
displacement and stress vary according to Bessel functions 
rather than trigonometric functions as far as the radial 
coordinate is concerned. For large values of r , we have 

_ , - sin * + as.ii , , - *m& — oa£7 \*J) 

JJ.&}-* " .J,'jTJ-> — 



Thus, far from the origin the motion becomes periodic in y . 
Asbecomes very large, the straight crested behavior is the 
limit of circular crested waves. 

The secular equation (28) is the transcendental equation 
and contains complete information about the phase velocity, 
wave number and attenuation coefficient of the plate waves. 
In general, wave number and hence the phase velocity of the 
waves is complex quantity, therefore the waves are attenuated 
in space. 

In order to find real phase velocity and attenuation 
coefficient, we write 



(26) 
R and q are real. 



c-' = V - isf'Q 
so that £ = R + i Q , where v = — 

* ^ R 

Here y is the propagation speed and Q the attenuation 
coefficient of the waves. Upon using relation (26) in (23) the 
values of y and Q for different modes have been obtained 
numerically. 
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VI. SPECIFIC LOSS 

The ratio of energy dissipated a W in taking a specimen 
through a cycle, to the elastic strain energy (w ) stored in 
the specimen when strain is maximum is called specific loss. 
It is the most direct method of defining the internal friction 
for material. According to Kaliski [2] for a sinusoidal plane 

AW 
wave of small amplitude, the specific loss 



II' 



equals 



4;r times the absolute value of imaginary part of E, to the 
real part of t, , i.e. 
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VII. UNCOUPLED THERMOELASTICITY 

In case of uncoupled thermoelasticity (UCT), the 
thermomechanical-coupling constant vanishes, i.e. e = 0, 
hence (15) reduces to 



m* + B m 2 + C =0 and m 



(27) 



where B 



$ 2 {p-Jc 2 ) r- r(c 2 - C 2 Xl-c 2 ) 
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Consequently the secular equation (23) in this case 
reduces to 
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where P q9 ,f q9 ,q = \5,l, P lq ,f lq , q =5,7, P ql ,f ql ,q =5 are 
given by (41) and (31.1). 
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m q , q = 1,5 are the roots of (28) 

For inviscid liquid (p. L = o ) , at uniform temperature we 
get the secular equation same as obtained and discussed by 
Sharma and Pathania [6]. 

VTA. NUMERICAL RESULTS AND DISCUSSION 

In order to illustrate the theoretical results obtained in 
the previous section, some numerical results are presented. 
The materials and liquid chosen for this purpose is zinc and 
water. From Sharma and Pathania [7] and Sharma and Sharma 
[9], we have taken physical data. 

/j=S.536xl0^fl* i , ^=l#lXy\# l Niri\ c l2 =1.650x10 ''Nm 2 

c :i =im7xl(}\Vm\ Cj^lS&lxld'xVflT, c M =1 ^tOxl^A^ 
fi =704xltf Km 1 As.g T fi =690xl0 i A r flB i deg„s=l29xl<f 2 

^ =0.69(kltf Wrt de|, K, =0.690dtfffirf deg r d =L5Sxltf L 5 H 

C m = ADJyzWJKgm- y fes ! d = 1,. h = 0.25 H = 01 

c L =1.5k10 j arfj ,p L =1000 Kgm* r, = 2$$ K 
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The variations of phase velocity of both symmetric and 
asymmetric modes of wave propagation for inviscid (ideal) 

liquid(p L ±0.0,/u L =0.0), low viscous liquid 
(p L * 0.0,1.1 L = 0.1) and high viscous liquid 
(p L ^0.0,^ =1.0) are plotted in Fig. 1 with respect to the 
wave number for angle of inclination 45 ° . 

It is observed that the phase velocity of lowest (acoustic) 
asymmetric mode (AJ is observed to increase from zero value 
at vanishing wave number to become closer to Rayleigh wave 
velocity at higher wave numbers in direction of wave 
propagation. 

The velocity of the acoustic symmetric mode (S ) become 
dispersion less i.e. remains constant with variation in wave 
number. The phase velocity of higher (optical) modes of wave 
propagation, symmetric and asymmetric; attain quite large 
values at vanishing wave number, which sharply slashes 
down to become asymptotically closer to the shear wave 
speed The velocity of higher modes is observed to develop 
at a rate, which is approximately n-time, the magnitude of the 
velocity of first mode (n=l). 

The surface wave propagates with attenuation due to the 
radiation of energy into the medium. This radiated energy 
will be reflected back to the center of the plate by the lower 
and upper surfaces. In fact, the multiple reflections between 
the upper and lower surfaces of the plate form caustics at 
one of the free surface and a strong stress concentration 
arises due to which wave field becomes unbounded in the 

limit d — > 00 .The unbounded displacement field is 
characterized by the singularities of circular tangent functions. 
The variation of attenuation coefficient with wave number 
for symmetric and asymmetric acoustic modes is represented 
in Fig. 2. 

The attenuation coefficient for acoustic mode has almost 
negligible variation with wave number. For higher mode of 
wave propagation, the magnitude of attenuation coefficient 
increases monotonically to attain a maximum value and then 
slashes down sharply to zero with the increase in wave 
number For first mode of wave propagation for high viscous 
the magnitude of attenuation coefficient shoots up then 
slashes down sharply to zero with the increase in wave 
number. 

Fig. 3 shows the variations of specific loss with wave 
number for symmetric and asymmetric modes for angle of 
inclination. The variations of specific loss factor of energy 
dissipation of symmetric and asymmetric acoustic modes are 
plotted with respect to wave number. 

The trend and behavior of symmetric and asymmetric 
profiles is noticed to be dispersive. The profile of asymmetric 
mode attains local maxima at R = 1 , 4 and minima near about R 
= 3 before it becomes smooth and stable for R>7. Fig. 4 
presents variations of specific loss factor of energy 
dissipation with respect to liquid loading temperature. 
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The specific loss factor profile of asymmetric acoustic mode 
remains non-dispersive and varies linearly with liquid loading 
temperature; the specific loss factor of asymmetric mode for 
(n=0) has small magnitude as compare to the symmetric 
mode. Further for other asymmetric mode the magnitude is 
greater than the symmetric mode. The asymmetric mode for 
(n=l), the magnitude is greater than (n=2). The profiles of 
specific loss factor of symmetric and asymmetric modes follow 
linear and non-dispersive trends at all values of the loading 
temperature. The magnitude of specific loss factor of 
asymmetric mode is noticed to be quite high as compared to 
that for all other considered cases. 




Figure 1 . Variation of phase velocity with wave number 




Figure 2. Variation of attenuation coefficient with wave number 
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Figure 3. Variation of specific loss with wave number 
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Figure 4. Variation of specific loss with liquid temperature 

CONCLUSIONS 

The effect of the viscous loading on the phase velocity, 
attenuation coefficient due to viscosity have been computed 
and discussed. It is observed that the inviscid liquid loading 
does not contribute much to system. However, increasing 
viscosity of liquid loading magnifies the value of phase 
velocity of both symmetric and asymmetric modes. The 
profiles of attenuation coefficient and specific loss factor of 
acoustic modes are noticed to be highly dispersive. 
Significant effect of liquid temperature has been observed 
on specific loss factor of energy dissipation profiles in the 
considered material plate. 
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